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Abstract 

We calculate tunneling conductance of normal metal and superconducting topological insulator Cu^I^Ses 
junctions for the possible pairing symmetries. In the presence of gapless Andreev bound states (ABSs), the 
tunneling conductance shows a zero-bias peak even for the three-dimensional full-gap superconducting state. 
This zero-bias conductance peak stems from the enhancement of the surface density of states induced by 
the surface-state transition in momentum space. 
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1. Introduction 

Topological superconductors (TSCs) are a new state of matter which is characterized by the nonzero 
topological numbers of the bulk wave functions 0,0, S|- F rom the bulk-boundary correspondence, there exist 
topologically protected gapless Andreev bound states (ABSs). In particular, the superconductivity infers 
that the gapless ABSs are their own antiparticles, thus Majorana fermions. The realization of Majorana 
fermions in condensed matter physics is of particular interest because of their novelty as well as the possible 
application for quantum computing. The recent discovered superconductor Cu a: Bi2Se3 [j] is a candidate of 
the TSC. The parent material Bi2Se3 is a topological insulator with a topologically protected gapless Dirac 
fermion on its surface. With doping Cu, however, the superconductivity appears. From the Fermi surface 
structure of the material, it was predicted that the superconducting topological insulator (STI) Cu x Bi2Se3 
must be a TSC (gj if time- reversal invariant odd-parity superconductivity is realized 0,01 • Possible ABSs 
specific to this material have been predicted @, [jj [l(J ■ Indeed, point-contact spectroscopy experiments 
[ill (HI have reported a zero-bias conduction peak (ZBCP) which signifies gapless ABSs [l3, 1J, HI- On 



the other hand, the superconducting proximity effect on topological insulators has been also studied [lQ, 
17, ID, The purpose of this work is to present a theory of the tunneling conductance for the STI. For 



one of the possible pairings, the ABS has a linear energy dispersion. In this case, SDOS does not always 
express the tunneling conductance in normal metal (N)/superconductor junction 2^, mm mils, and 



the resulting tunneling conductance depends on the transparency at the interface. Moreover, the BW state, 
which is a three-dimensional full-gap TSC and supports dispersive gapless ABSs, never shows a ZBCP 



25| . The relation between SDOS in three-dimensional full-gap superconductors and ZBCPs is non-trivial. 
Thus, in order to pursue the origin of the ZBCP in the STI, one needs to calculate directly the tunneling 
conductance of the N/STI junction. In this work, we calculate the tunneling conductance between N/STI 
junctions for all possible pairing symmetries. The STI with dispersive gapless ABSs can show a ZBCP, 
due to the enhancement of the SDOS induced by the surface-state transition in momentum space. This 
transition occurs since the parent material (a topological insulator) has a surface Dirac fermion which still 
remains after Cu doping. 
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Table 1: Momentum- independent paring potentials Aj in STIs with point group symmetry [Eg]. The irreducible represen- 
tation, parity eigenvalue, and structure of the energy gap are indicated. 
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2. Model 

A model Hamiltonian of a superconducting topological insulator is given by 

H (k) = -fj, + m(k)a x + v z k z a y + va z (k x s y - k y s x ), (1) 
H(k)=H Q {k) Tz +A t T x , (2) 
m(fe) = mo + mik z + to 2 (k x + k y ) , (3) 

where the basis is set to be (<V|-(fc), c CT j_(fc), — .(—fc), c^(— fe)) T , cr = ± and s =|, | denote the orbital and 
spin, (Tj, and are Pauli matrices in the orbital, spin, and Nambu spaces, respectively. This system can 
have four types of momentum-independent pairing potential A, [5[ , which are summarized in Table [TJ 

The energy gap for each pair potential is shown in Fig. 1. For Ai a , an isotropic full-gap superconductor 
is realized. On the other hand, for Aib and A 2 , the energy gap fully opens but is anisotropic. For A 3 and 
A4, point nodes appear on the pole (fc z -axis) and on the equator (k x or k y -axis), respectively. 



3. Surface-state transition in the momentum space 

On the surface perpendicular to the z-axis, STIs with A 2 and with A 4 support gapless ABSs @, i, i, [loj] . 
Figure [5] shows energy spectra of the gapless ABSs of STIs with A 2 and A4 for different values of m\. There 
are four types of ABSs, i.e., (a): cone type for A 2 , larger \x and toi, whose energy spectrum is monotonic 
function of k, (b): caldera type for A 2 , smaller \x and toi, where the energy spectrum has zeros both at 
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Figure 2: Energy spectra of the gapless ABSs. Only the half is shown. The whole energy spectrum is given by mirroring with 
respect to E = 0, preserving particle- hole symmetry. Pairing potentials are set to A 2 for (a) and (b), and to A4 for (c) and 
(d). For (a) and (c) [(b) and (d)] we choose mi = 20.18eVA 2 (mi = 5.65eVA 2 ), where cone and ridge (caldera and valley) 
type energy spectra are realized, respectively, (e): The energy spectra as a function of kit/hp with feii = (fc 2 + fc 2 ) 1 / 2 for A2 
near the transition point (mi ~ mi c = mot> 2 //^ 2 ). k-p denotes the Fermi momentum for fc z = 0. 



k = and at finite k, (c): ridge type for A4, larger /1 and mi, where zero energy states appear between the 
two point nodes (on the fc y -axis), (d): valley type for A4, smaller /i and m±, where the energy remains zero 
on the fc y -axis but becomes a monotonic function of k x and zero-energy states appear at a finite momentum. 

Differently from the other topological superconductors, these STIs can have twisted shape of energy 
spectra of gapless ABSs @,@> E3] for smaller /1 and mi, as shown in Figs, ^h) and[2Jd). Zero-energy states 
at finite momenta stem from surface Dirac fermions in the normal state. In other words, pairing potentials 
of A2 and A4 cannot open superconducting gap for the surface Dirac fermions, protected by the mirror 
symmetry (jj [l(|, as explained below. On the other hand, for much larger /x and mi where the Fermi level is 
much higher than the energy of surface Dirac fermions, the resulting energy spectrum of the gapless ABSs 
are not twisted since surface Dirac fermions are irrelevant in this energy regime. Consequently, the transition 
from the twisted one to the conventional one occurs at the intermediate values of /i and mi as mi/i 2 = mov z 
[ioj . Then, at this transition point, the energy spectrum of the ABS is proportional to fc 3 [Fig. H[e)] and 
the corresponding SDOS diverges. In the vicinity of the transition point, the tunneling conductance yields 
a ZBCP due to the divergent SDOS at the transition point, as we will see in the next section. 

Here, we explain symmetry-protection of the zero-energy ABS at a finite momentum for A2, in terms of 
the mirror symmetry |9l. llOj. The system has mirror symmetry as s x T z Ho(—k x , k y , k z )T z s x = Ho(k x , k y , k z ). 
Then, s x t z = ±1 becomes a good quantum number for k x = 0. The particle [Ho(k)] and hole [— Ho(k)] 
components of a surface Dirac fermion for A = have the different mirror eigenvalues. On the other hand, 
A2 = Aa y s z is even under the mirror reflection as s x t z /S.2T z s x — A2. As a result, due to mismatch of 
the mirror eigenvalues, A2 cannot open the superconducting gap of the surface Dirac fermions. Namely, 
when the system has surface Dirac fermions at the Fermi level, they still remain gapless and turns into 
zero-energy ABSs. In this picture, the zero-energy ABSs are located at the Fermi momentum of the surface 
Dirac fermion &df m the normal state. As one increases fi and mi, &df decreases. For certain values of 
/j, and mi with &df = ^f, the surface Dirac fermions are gone from the Fermi level. However, the ABSs 
remains gapless since the ABSs with the different mirror eigenvalues are not hybridized. Further increasing 
/j, and mi, the momentum of the zero-energy ABSs decreases to 0, then the shape of the energy spectrum 
changes from twisted to non-twisted, i.e., the surface-state transition occurs. Note that the zero-energy 
ABSs protected by the mirror symmetry on the fcj,-axis are robust even if the higher order terms of k, e.g., 
the warping term 26], is taken into account. 
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Figure 3: Tunneling conductance as a function of bias voltage V. Gn/Go denotes the normalized conductance of the normal 
state (A = 0), which corresponds to the transmissivity. The parameters of STI are the same as in Fig. [2] and fJ.fi/ ^ = 0.6, 60, 
and 1200 for Gn/Go = 0.88,0.008, and 0.0001, respectively. Effective mass of electron in the normal metal is taken to be 
m N = l/m 2 . 



4. Tunneling conductance 

Now, we discuss the tunneling conductance of N/STI junctions. Hamiltonian of the normal metal is set 
to Hjsi{k) = (k% + ky + k 2 )/{2m^) — /in- We calculate the tunneling conductance by using BTK formula [Io| 
for various /zn, which corresponds to transmissivity of the junction. When fi^/fi *~ 0.6, the transmissivity 
of the junction is high Gn/Go ~ 1 since the magnitudes of the Fermi momenta of the normal metal and of 
the STI are the same. Here, conductance of a perfect transmission is given by Go = k"^ F L 2 / (4ir) x e 2 /h, 
where ^nf is Fermi momentum of the normal metal and L 2 denotes cross section of the system. As one 
increases fiist, the magnitude of Fermi momentum of the normal metal increases and then the transmissivity 
decreases. 

Figures [3] and H] show the tunneling conductance G normalized by that of the normal (A = 0) state Gn- 
The line shapes of the tunneling conductances for Ai a and Aib are the same as that of the conventional 
s-wave superconductor, as shown in Figs. EUa) and [21b). For A3, the system has point nodes on the fc z -axis 
[Fig. [TJd)] and has no gapless ABS on the surface perpendicular to z-axis. The corresponding tunneling 
conductance is given by G oc V 2 with V being bias voltage [Fig. EHc)]. 

On the other hand, for A2 with a conventional cone type ABS [Fig. Ufa)], a ZBCP is realized for the 
intermediate transmissivity (Gn/Go = 0.008). It is noted that a ZBCP never appears for the BW state 
[25| . which is a well-known three-dimensional TSC. In the STI, as discussed in the previous section, the 
SDOS diverges at the zero-energy for the surface-state transition point. Therefore, in the vicinity of this 
transition point, the magnitude of SDOS near the zero-energy becomes large. This results in a ZBCP 
in the intermediate transparent case. In the low transparent limit (Gn/Go = 0.0001), the line shape of 
the tunneling conductance completely converges to that of the SDOS. Thus the line shape of tunneling 
conductance for Gn/Go = 0.0001 forms a double-peak and a zero-bias dip [Fig. Ufa)], which is the same as 
that of the BW state. For A2 with the caldera-type ABS, the corresponding SDOS near the zero energy 
becomes larger hence the height of the resulting ZBCP becomes larger [Fig. Etb)]. In addition, the ZBCP 
survives in a lower transmissivity (Gn/Go = 0.0001). 

Finally, for A4, the resulting tunneling conductance shows a ZBCP for any transmissivity [Figs. HIc) and 
Hid)], since zero-energy ABSs appear between the two point nodes, i.e., a fiat band on the k y -axis appears. 
For smaller mi and /i, the gapless ABSs changes from the ridge to caldera types. Similarly to that of caldera 
type ABS of A2 [Fig. HIb)], the resulting tunneling conductance becomes prominent [Fig. Hid)]. 

5. Summary 

In this work, we have clarified tunneling conductance of normal metal / STI junctions for possible pair 
potentials. Only STIs with A2 in the intermediate transmissivity and A4 show a ZBCP. The presence of 
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Figure 4: Tunneling conductance and the corresponding gapless ABSs of the STIs. The parameter are the same as in Figs. [2] 
and [3] 

ZBCP is possible in the STI with A2 although the energy spectrum of gapless ABS is dispersive. A ZBCP 
in a STI with A2 originates from the surface-state transition, which comes from surface Dirac fermions in 
the normal state, i.e., topological insulator. This is not the case in the BW state. Our obtained results serve 
as a guide to explore three-dimensional topological superconductor hosting Majorana fermions. 
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